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The effect of ambient disorders and sequence heterogeneities on the reptation of a long polymer is 
studied with the aid of a disordered tube model. The dynamics of a random heteropolymer is found 
to be much slower than that of a homopolymer due to collective pinning effects. The asymptotic 
properties belong to the universality class of a directed path in (f +1 )-dimensional random media. 



PACS: 83.20.Fk, 83.fO.Nn, 05.60.+j, 36.20.Ey 

The stochastic motion of a polymer chain entangled in 
a disordered medium such as a gel is of much scientific 
and technological interest. A convenient starting point 
for understanding the polymer dynamics is de Gennes' 
reptation model, which describes the polymer's wormlike 
motion along a fictitious tube threaded through an array 
of fixed obstacles |jf|-§| . In the original reptation model, 
the polymer was assumed to be homogeneous, and effects 
due to randomness in the medium (e.g., spatial varia- 
tions in the pore size of a gel) were neglected. A later 
argument by Harris Q showed that the static configura- 
tion of a self-avoiding polymer is in fact not affected by 
the randomness, a result supported by numerical simula- 
tions However, extensive numerical [^|-^| and exper- 
imental |],[l^] studies found the dynamics of polymers in 
random media to be much slower than the classical repta- 
tion dynamics Jl] || • These results are understood quali- 
tatively within an "entropic trapping" framework H[n] . 

It should be noted that although some of the best 
known applications of reptation theory are concerned 
with the behavior of biopolymers which arc inherently 
heterogeneous, there has been little theoretical work 
on the reptation of heteropolymers beyond a restricted 
model analyzed by de Gennes many years ago |l2| . In 
this article, we address the combined effect of the het- 
erogeneity of the polymer and the randomness of the 
media. Polymer motion is separated into two compo- 
nents: reptation within a tube and diffusive motion of 
the center-of-mass. For a homopolymer, the disorder 
does not affect the reptational motion, but the c-o-m 
motion is drastically slowed down due to entropic trap- 
ping. For a heteropolymer, the reptational motion is also 
drastically affected: The well-known algebraic reptation 
time becomes exponentially long, making the dynamics 
of the heteropolymer even much slower. Our results are 
obtained with the aid of a disordered tube model, which 
we find to be analogous to the well-known problem of a 
directed path in random media [jp| . 

Consider first a long, flexible homopolymer entangled 
in a covalently-crosslinked gel matrix, characterized by 
a typical pore size a. We shall take a > £ ^S> b, £ 
being the persistence length of the polymer and b be- 
ing the monomer size. (E.g., double-stranded DNA in 



agarose gel has a ~ 5000,4, £ » 500A and b = 3AA @; 
pore sizes are typically much smaller in polyacrylamidc 
gels . ) Since the equilibrium configuration of the poly- 
mer is that of a self-avoiding walk, described by the size 
exponent v > 1/2 even in the presence of disorders jfp|§, 
we shall describe the dynamics of the polymer by the 
reptation tube model ||. In this model, the polymer 
is confined to a tube of diameter ~ a |], and mod- 
eled as "beads" linearly connected by (entropic) springs 
of spring constant K = §/cbT/£ 2 . Configurational en- 
tropy associated with the exponential number of (self- 
avoiding) tube trajectories results in an effective tensile 
force acting on the two ends of the bead-spring chain 
as in the disordered-free case [|. This effect is incor- 
porated into the model by letting the springs acquire a 
finite equilibrium length of the order a |^]. Each link of 
this bead-spring chain therefore represents an elementary 
"blob" of Mbiob ~ fl2 / (bQ monomers whose physical size 
is ~ 0(a). Confinement of the polymer to the tube costs 
a free energy (per blob) of the order V = fc^T • (£/a) 2 [||. 
Since the diameter of the reptation tube embedded in 
the disordered gel matrix is non-uniform along its length, 
we describe the variable confinement entropy by a "ran- 
dom potential" V(s), where s denotes the curvilinear 
coordinate along the tube. For a typical self-avoiding 
tube trajectory, V may be modeled by a short-range cor- 
related random variable, characterized by the variance 
6V(s)6V(s') = A v S a (s - s'), where 5V(s) = V(s) - V, 
6 a (s) is an exponentially damped function of range a, and 
the overline denotes average over the ensemble of tubes. 
Assuming that the fractional variance in tube diameter 
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variation is O(l), we have Ay ~ V . 

A key advantage of the original tube model Q is that 
the nonlocal cxcluded-volume interaction between the 
beads can be neglected. This results from the entropy- 
generated stretching forces, which makes the chain ex- 
tended in the tube, i.e., the contour length of the chain 
being linearly proportional to the number of beads. Con- 
sequently, the large scale dynamics of the polymer is ob- 
tained simply as the Id Rouse dynamics of the beads. 
Let us assume that the nonlocal interaction between the 
beads (confined in the self-avoiding tube) can be ne- 
glected also for the disordered case. Then the statisti- 
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cal mechanics of a long polymer with degree of polymer- 
ization M 3> Mbiob is given by a bead-spring chain of 
N ~ M/Mbiob 3* 1 links, with the "Hamiltonian" 

N (K ~) 
^o = E yK+i-s„-a) 2 + y( Sn )L (1) 

n=0 ^ J 

where s„ denotes the coordinate of the n th bead along 
the tube. 

The legitimacy of neglecting the nonlocal bead- 
bead interaction certainly deserves scrutiny: As in the 
disorder-free case, the approximation is justified provided 
the chain is in an extended conformation. On the other 
hand, a Gaussian chain in a tube can in principle col- 
lapse |1| to a region where SV is large and negative, 
corresponding to a section of the tube with wider open- 
ings. If the chain described by (|l|) collapses, then the tube 
model would not be self-consistent. What opposes the 
collapse of the Gaussian chain is again the tensile stretch- 
ing forces generated by the configurational entropy of the 
self-avoiding tube. This physics is reflected in the model 
(|l|) as a competition between the "elastic energy" cost 
against chain collapse, of the order ja 2 per spring, and 
the disorder energy gain, of the order \/ Ay • The out- 
come of this competition has been obtained recently in 
the context of "non-Hermitian quantum mechanics" |L6| : 
A discontinuous phase transition between the collapsed 
and the extended state is found at a certain critical point 
given by the condition -|-a 2 ~ V Ay- Using the expres- 
sions given above for K and Ay for the polymer prob- 
lem at hand, we find that the extended state is always 
preferred under the presumed condition a > £. Thus, the 
disordered tube model (|l|) is justified self-consistently. It 
will be used from here on to generate the large scale dy- 
namics of self-avoiding polymers in random media. 

Classical reptation dynamics is characterized by the 
reptation time tr.q ~ N 3 , which is the time it takes for 
the polymer to reptate from a given tube to a completely 
different tube in the adjacent neighborhood. For time 
scales much exceeding tr^, the polymer behaves as a 
point particle undergoing Brownian motion. This leads 
to a diffusion coefficient D (N) ~ N 2v /r Rfi (N). In the 
presence of randomness, the reptation time tr is modified 
(see below); but more significantly, the large scale mo- 
tion becomes one of thermally-activated barrier hopping, 
since adjacent tubes may be narrower and hence higher 
in energy The typical barrier height Eb(N) can be 
taken as the free energy variation of a polymer confined 
to cells of size ~ N v . This variation is bounded to be 
within the order of iiV 1 / 2 and is found numerically |l7]] 
to scale as TV", with u> w 0.15 in 3d and u> « 0.28 in 
2d. Thus the overall diffusion coefficient is reduced dras- 
tically to D ~ [N 2l/ /T R (N)]e- Eb{N}/kBT for large N. 

To estimate the reptation time tr(N) itself for the 
random system, it is necessary to separate the sliding 
motion of a polymer along its prescribed tube from the 



thermally-activated barrier hopping. This can be accom- 
plished within the disordered tube model ([!]) by artifi- 
cially imposing a periodic boundary condition on the ran- 
dom potential V, i.e., V(s) = V(s + Na). With this 
simplification, the extended state of the chain can be ob- 
tained straightforwardly, by observing a discrete trans- 
lational symmetry s„ — > s rl +i which the model pos- 
sesses even in the presence of disorder: The polymer can 
move along the tube by the propagation of a longitudinal 
"defect" or "kink" (e.g., a configuration with s n = 
much like what was proposed in the original reptation 
model [0. Thus, the random potential is irrelevant in 
the extended state Jl6| , and the scaling properties in the 
sliding regime are the same as those in the disorder-free 
case. For example, tr ~ tra, accompanied by a contour 
length fluctuation SL {N) w {{k B T / K)^ 1 ' 2 . 

The above analysis clarifies the role of entropic trap- 
ping for a homopolymer: The very slow motion of a long 
homopolymer in random media is dominated by the ex- 
ponentially long waiting time ~ tr } q e N / fcsT needed to 
overcome large spatial variations in the mean confine- 
ment potential. It is not a result of local "bottlenecks" 
which impedes the reptational motion along the tube. 
To study this behavior experimentally or numerically, it 
is necessary to have sufficiently long polymers such that 
each polymer threads through many local bottlenecks. 

We next consider a heteropolymer which may carry 
side groups of vastly different sizes, or may have different 
degrees of partial charges, hydrophobicities, etc. Follow- 
ing de Gennes jl2| , we model the coarse-grained effect of 
the side groups by assigning a "charge" q n 6 {±<2} to 
each bead n. Then the interaction energy V in Eq. (g) 
becomes explicitly n-dependent, with the form 

V(s n , n) = V (s n ) + q n Vi(s n ), (2) 

where Vo(s) ± QVi(s) is respectively the free energy cost 
of placing a bead of ± charge at the coordinate s along 
the tube. We will take the random potentials Vq and 
V\ to be short-range correlated, with variances of the 
order Ay. To model a heteropolymer, we choose the 
charges {q n } randomly, with probability p for q = +Q 
and 1 — p for q = —Q, such that [q] = (2p — 1) Q, and 
[Sq m Sq„] = A q 5 m , n , A q = 4p(l-p)Q 2 . (Here, [...] 
denotes average over the ensemble of random sequences.) 
The magnitude of Q depends on the coarse-graining scale 
a and variations of the monomer-tube interaction. Let 
the latter be characterized by a fractional variance o~, 
then we have Q 2 = crMbiob- 

As in the homopolymer case, a Gaussian chain de- 
scribed by the Hamiltonian ([!]) with the interaction (j^) 
may be in cither the collapsed or the extended state de- 
pending on the magnitude of the elastic energy ^a 2 . 
However for a RHP, different monomers tend to prefer 
different sections of the tube. Hence the extended state 
is even more preferred, and we are again justified to ne- 
glect the non-local bead-bead interaction. 
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Let us consider now the reptational motion of an ex- 
tended RHP, again by imposing periodic boundary con- 
dition on the V's. Heterogeneity in sequence composition 
has profound effects on the dynamics, as the random q„'s 
remove the translational symmetry s n — ► s n +i described 
above for the homopolymer. Thus, larger beads prefer 
to occupy segments of the tube with wider cavities, etc. 
De Gennes demonstrated the existence of such a het- 
cropolymeric effect by considering a singular limit of the 
RHP model [E| , where the longitudinal elasticity of the 
polymer is suppressed by taking K — > oo. Below, we 
will show that a significant heteropolymer effect can ex- 
ist generically for an elastic RHP. 

To proceed, we introduce a longitudinal displacement 
field, u n — s n — na, and write the Hamiltonian as 



H = \ y K+i - u n ) 2 + W{u n , n) \ 

n=l ^ > 



(3) 



where W(u n ,n) = 5q n ■ 5V\(u n + na,n). [There are 
actually two additional terms in TL: q n Vi and 

S„ (Vb(sn) + [?]^(«ti))- The first term is simply an over- 
all energy shift and does not affect the motion of the 
polymer along the tube. The second term is just like 
the random potential of the homopolymer problem p); 
it is irrelevant as discussed above.] The Hamiltonian (g) 
then describes a fictitious "directed path" with "trans- 
verse" coordinates {«„}, embedded in a 1 + 1 dimen- 
sional random medium W(u,n) [fL8| . The randomness 
results from a combination of sequence heterogeneiity 
and medium disorder. To elucidate the properties of this 
system, we first perform a naive perturbative analysis 
for small randomness A = A 9 Ay. We find the ran- 
domness to have a negligible effect on the classical rep- 
tation results for chains below a crossover scale N x w 
(kBT) 5 /(KA 2 a 2 ). The perturbative analysis fails for 
long chains with iV > N x , indicating that the asymptotic 
reptation properties are qualitatively affected by arbitrar- 
ily weak sequence heterogeneities. Note the dependence 
of the polymer's crossover length M x = N x Mbiob on the 
tube size a: Using expressions for K and A given above 
for the simple RHP model, we find M x ~ a 4 b/a 2 ^ 5 , 
which increases quickly with increasing a/£ as suggested 
in Ref. E^], but is nevertheless accessible for sufficiently 
heterogeneous polymers in typical gels. 

To obtain the asymptotic behavior of the polymer with 
M ^> M x , regime, let us examine more closely the effec- 
tive random potential W(u,n): The statistics of W is 
easily obtained in terms of the statistics of q n and Vi(s), 
with JW] = and 



[W(u, n)W(u', n')] =AS a (u- u')8„ 



(4) 



The correlator (^) indicates that an effective two- 
dimensional random "point" potential is generated, even 
though W itself, being a product of two one- dimensional 



random variables, must contain long-range correlations. 
The latter is manifested in the higher moments of W . For 
instance, W 2 {u,n) contains a term [(<5g) 2 ] SV 2 (u n + na) 
which is correlated along the direction u n = —na. But 
such terms are just like the random potential V(s) of 
the homopolymer case and are irrelevant. Thus, we con- 
jecture that the system ^) belongs to the universality 
class of a directed path in (1+1) -dimensional uncorre- 
cted Gaussian random potential. 

This universality class is well known |l3|]. In the 
asymptotic regime, the polymer exhibits glassy dynam- 
ics which is characterized by two exact scaling laws: 
The sample-to-sample free energy variation scales as 
SF(N) w k B T ■ (iV/iVx) 1 / 3 , and the contour length fluc- 
tuation scales as SL(N) ~ s x ■ (N/N x ) 2 / 3 , where s x — 
SL (N X ) = (k B T) 3 /(KAa). The conjectured equiva- 
lence between RHP reptation and the directed path prob- 
lem were tested, by numerically computing SF(N) and 
SL(N) from the model (0) and the interaction (Q), treat- 
ing K , A as arbitrary parameters. We used the transfer 
matrix method |L3) , and examined particularly the "zero- 
temperature" limit of the model, whose small crossover 
length i¥ x enables us to access the asymptotic scaling 
regime. The zero-temperature problem is an optimiza- 
tion problem. It is necessary to place the bead-spring 
chain on a discrete lattice. To simplify the numerics, we 
restrict the displacement of beads such that s n+ i — s n 
can only take on the values {a, a ± 1} at each step n of 
the transfer matrix. A random potential V±(s) of zero 
mean and unit variance is assigned to each lattice point, 
and a charge q n is assigned independently to each bead. 
The numerical results for the disorder-averaged contour 

1/2 

length fluctuation SL — [{u* N ) 2 } of the optimal con- 



figuration u* (with the n = end fixed at the origin) 
and the sample-to-sample variation SE*(N) of the to- 
tal energy of the optimal configuration E*(N) [jl9| are 
shown in Fig. 1(a). It is seen that fluctuations rapidly 
approach those expected of the directed path, 5L ~ TV 2 / 3 
and SE* ~ N 1 ^ 3 , suggesting that the RHP is indeed in 
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FIG. 1. Heteropolymers confined in a tube: (a) Scaling of 
the contour length fluctuation 8L, and the sample-to-sample 
minimal energy variation 8E* , for bead-spring chains of 
N = 4000, a = 4, K = 1, and [q] = 0.1, averaged over 
1500 samples, (b) Auto-correlation of the energy landscape 
E(s,N); the straight line has slope 1. 
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the same universality class as the 1+1 dimensional di- 
rected path. Similar behaviors have recently been found 
in a number of related studies |^,|lj . 

Another useful quantity to examine is variation in the 
polymer's "energy landscape" , which we obtain by fixing 
the n = end of the polymer to an arbitrary coordi- 
nate s, and then computing the energy E(s, N) of the 
optimal configuration (of the constrained polymer) for 
each s. The landscape is characterized by the correlation 
function C E (s -s',N)= \(E(s, N) — E(s', AO) 2 ] , which 
is expected [O] to have the scaling form Ce(s,N) = 
N 2 / 3 g(\s\/N 2 /•*), with the scaling function g(x % 1) <~ x 
and g(x > 1) ~ const. This scaling form is verified by 
our numerics as shown in Fig. 1(b). 

The zero-temperature behavior described here corre- 
sponds to the asymptotic regime of the RHP, beyond 
the crossover scale -/V x . For the finite temperature prob- 
lem at hand, the RHP dynamics along the tube is gov- 
erned by variations in the free-energy landscape F(s, N), 
whose form is obtained easily from E(s,N). It is de- 
scribed by the correlation function C F (s — s',N) — 
{k B T) 2 \s- s'l/sx, for s x < \s-s'\ < SL{N), satu- 
rating to C F ~ (SF) 2 w (k B T) 2 (N/N x ) 2 / 3 for larger 
displacements \s — s'\. 

Dynamics of the chain in such a rough free energy 
landscape cannot proceed by the propagation of a few 
isolated kinks. It requires instead collective movement 
of large pieces of the chain such that the chain always 
remains in local optimal configurations |2^| . The re- 
sulting "collective creep" dynamics can be modeled by 
the motion of one of the chain ends in the Id Brow- 
nian random potential F(s,N) |25(] . The latter is a 
well-studied classic problem p£|. The time it takes for 
the particle to move by a distance s > s x along the 
tube is T Rfi (N x ) ■ exp [C F (s)/(k B T) 2 ] ~ e s / s *. Since 
C F {s)/(k B T) 2 saturates at (N/N x ) 2 ? 3 for 6L < s < L, 
the reptation time is 



tr(N) « r fl „ (^x)exp[(iV/iV x ) 2 / 3 ], 



(5) 



which greatly exceeds the classical reptation time 
tr,o(N) if the polymer is sufficiently long. Note that tr 
also exceeds the waiting time ~ e NUJ / k BT neec } ec [ f or 
activated barrier hopping by the homopolymer. Hence 
we conclude that the dynamics of a long random het- 
eropolymer is much slower than that of a homopolymer. 

The extraordinary slow dynamics of the RHP is a col- 
lective phenomenon resulting from a global "resonance" 
of the randomness in the polymer composition with the 
randomness in the tube structure; it is not a local effect, 
say, an increase of the microscopic friction. This can be 
tested by comparing an RHP with a ABAB copolymer 
consisted of the same monomers. Upon coarse-graining 
by a finite scale, the copolymer becomes equivalent to a 
homopolymer, while heterogeneities of an RHP cannot 
be coarse-grained away. Thus we expect a long RHP to 



exhibit much slower dynamics than a periodic copolymer 
of the same length and monomer content. 
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